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We study thermodynamic properties of the QCD vacuum in a magnetic field below chiral phase 
transition. The hadronic phase free energy in a constant homogeneous magnetic field is calculated 
in the framework of the chiral perturbation theory at non-zero pionic mass. It is demonstrated 
that the order parameter of the chiral phase transition remains constant provided temperature 
and magnetic field strength are related through obtained equation (the phenomenon of "quark 
condensate freezing”). 


I. INTRODUCTION 

The investigation of the vacuum state behavior under the influence of the various external factors is known to be one 
of the central problems in quantum field theory. In the realm of strong interactions (QCD) the main factors are the 
temperature and the baryon density. At temperatures below the chiral phase transition, T < T c , the dynamics of the 
system is characterized by confinement and spontaneous breaking of chiral symmetry (SBCS). At low temperatures, 
T < T c , the partition function of the system is dominated by the contribution of the lightest particles in the physical 
spectrum. In QCD this role is played by the n - meson which is the Goldstone excitation mode in chiral condensate. 
Therefore the low temperature physics (the hadron phase) enables an adequate description in terms of the effective 
chiral theory Ml- A very important problem is the behavior of the order parameter (the quark condensate ( qq )) 
with the increase of the temperature. In the ideal gas approximation the contribution of the thermal pions into 
the quark condensate is proportional to T 2 Si- In chiral perturbation theory (ChPT) the two - and three-loops 
contributions (~ T 4 and ~ T 6 correspondingly) into (qq) have been found in and II- At non-zero quark mass the 
analytic temperature dependence of the quark condensate was found in || in perfect agreement with the numerical 
calculations obtained at the three-loop level of the chiral perturbation theory with m q 0 j?]] . 

The situation with the gluon condensate ( G 2 ) = ((G“„) 2 )is very different. The gluon condensate is not an order 
parameter in phase transition and it does not lead to any spontaneous symmetry breaking (SSB). At the quantum 
level the trace anomaly leads to the breaking of the scale invariance and this in turn results in nonzero value of 
( G 2 ). However, this is not a SSB phenomenon and hence does not lead to the appearance of the Goldstone particle. 
The mass of the lowest excitation (dilaton) is directly connected to the gluon condensate, tod oc ((G 2 )) 1 / 4 . Thus 
the thermal excitations of glueballs are exponentially suppressed by the Boltzmann factor ~ exp {—m g i/T} and their 
contribution to the shift of the gluon condensate is small (A(G 2 )/(G 2 ) ~0.1 % at T =200 MeV) ||. Next we note 
that in the one-loop approximation of ChPT pions are described as a gas of massless noninteracting particles. Such 
a system is obviously scale-invariant and therefore does not contribute into the trace of the energy-momentum tensor 
and correspondingly into (G 2 ). It was shown in [[Io| that only in the three-loop approximation of ChPT does the gluon 
condensate become dependent on temperature in the chiral limit. At non-zero quark mass the analytic temperature 
dependence of the gluon condensate was found in [||. The effect of thermal excitations of massive hadrons on the 
properties of the quark and the gluon condensate was studied in [ |lT| within the generalized nonlinear cr-model invariant 
under conformal transformations. 

Another interesting problem is the study of the vacuum phase structure under the influence of the external magnetic 
field H. Quarks play an active role in shaping the QCD vacuum structure. Being dual carriers of both ’color’ and 
’electric’ charges they also respond to externally applied electromagnetic fields. The vacuum of strong interactions 
influences some QED processes has been discussed in Ref. Jl2|] . The behavior of the quark condensate in the presence 
of a magnetic field was studied in Nambu-Iona-Lasinio model earlier |h|. For QCD, the analogous investigation in the 
one-loop approximation of ChPT was done in []l4[ . It was found that the quark condensate grows with the increase of 
the magnetic field H in both cases. It implies that a naive analogy with superconductivity, where the order parameter 
vanishes at same critical field, is not valid here. The behavior of the gluon condensate (G 2 ) in the Abelian magnetic 
field is also a nontrivial effect. Gluons do not carry electric charge; nevertheless, virtual quarks produced by them 
interact with electromagnetic field and lead to the changes in the gluon condensate. This phenomenon was studied 
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in [^5|, jl(| based on the low-energy theorems in QCD. The vacuum energy density, the values of (G 2 ) and ( qq ) as 
functions of H have been found in the two-loop approximation of ChPT in |l6| . 

The low-energy theorems, playing an important role in the understanding of the vacuum state properties in quantum 
field theory, were discovered almost at the same time as quantum field methods appeared in particle physics (see, for 
example Low theorems jl7]]).In QCD, these theorems were obtained in the beginning of eighties [^8|. These theorems, 
being derived from the very general symmetrical considerations and not depending on the details of confinement 
mechanism, sometimes give information which is not easy to obtain in another way. Also, they can be used as 
” physically sensible” restrictions in the constructing of effective theories. An important step was made in Jl^|, where 
low-energy theorems for gluodynamics were generalized to finite temperature case. A relation between the trace 
anomaly and thermodynamic pressure in pure-glue QCD was obtained in |2(i|| by making use of the dimensional 
regularization in the framework of the renormalization group (RG) method. Also within the RG method, but by 
employing slightly different techniques an analogous relation was derived in the theory with quarks in |2l[| . 

The QCD phase structure in a magnetic field at finite temperature was investigated in ^2j. Within the ChPT 
framework the dependence of the quark condensate on T and H was studied and it was shown that the shift of 
the condensate is not the simple sum of the temperature (~ T 2 /F 2 ) and magnetic (~ H/(AttF^) 2 ) contributions. 
There appears an additional term which physically stems from the orbital diamagnetism of the charged pion gas. In 
the chiral limit the thermodynamic quantities, in particular the quark condensate, are characterized by additional 
dimensionless parameter A = y/H/T. However as it was shown in |^] the Gell-Mann-Oakes-Renner relation is not 
changed in a magnetic field at T ^ 0 (no additional terms appear) and thus the soft chiral symmetry breaking scenario 
remains the same. 

It should be noted, that ablian type chromomagnetic fields in QCD are essential for the nonperturbative dynamics 
of vacuum. The role of chromomagnetic fields in the thermal deconfining phase was investigated in |24[ ] . The influence 
of chromomagnetic vacuum condensate on the phenomenon of color superconductivity was studied in [p5| . The 
rearrangement of quark-antiquark vacuum in the strong gluonic fields was studied in |2fi| ] . 

In the present paper the vacuum free energy in magnetic field at finite temperature is calculated in the framework of 
ChPT. The general relations are established which allow to obtain the dependence of the quark and gluon condensates 
on T and H. A new phenomenon is displayed, namely the ’’freezing” of the chiral phase transition order parameter 
by the magnetic field when the temperature increases. The physical meaning of this fact is discussed. 


II. RENORMALIZATION GROUP PROPERTIES OF THE GLUON CONDENSATE ATT^O AND 

H + 0 

For non-zero quark mass (m q ^ 0) and in a magnetic field the scale invariance is broken already at the classical 
level. Therefore the pion thermal excitations would change, even in the ideal gas approximation, the value of the 
gluon condensate with increasing temperature and field 1 . To determine this dependence use will be made of the 
general renormalization and scale properties of the QCD partition function. 

The QCD Euclidean partition function with two quark flavors in external Abelian field A M has the following form 
(T = 1/P) 

Z = expi^-^ J dxi J d 3 xF 2 y | J [DB] [Dq] [Dq\ exp j - J dx 4 J c? 3 .t£ j . (1) 

Here the QCD Lagrangian in the background field is 

£ = ^2 ( G ^) 2 + 55 - iQqAf, - iyB“) + m 0q ]q, (2) 

q=u,d 

where Q q - is the matrix of the quark charges for the quarks q = ( u,d ), and for the simplicity the ghost terms have 
been omitted. 

The free energy density is given by the relation PVF (T, H, mo u , mod) = — In Z. Eq. (||) yields the following 
expression for the gluon condensate ((G 2 ) = ((G“j,) 2 )) 

„ gp 

(G 2 )(T, H , m 0u , mod) = 4 2 . (3) 

dip/go) 


1 


At zero quark mass and in the absence of magnetic field the gas of massless noninteracting pions is obviously scale-invariant and 
therefore does not contribute to the trace of the energy-momentum tensor and correspondingly to the gluon condensate ((Gj^) 2 ). 
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The system described by the partition function (fij) is characterized by the set of dimensionful parameters 
M,T, H,mo q (M) and dimensionless charge g^{M), where M is the ultraviolet cutoff. On the other hand one can 
consider the renormalized free energy Fr and by using the dimensional and renormalization-group properties of Fr 
recast (^]) into the form containing derivatives with respect to the physical parameter T, H and renormalized masses 
m. q . 

The phenomenon of dimensional transmutation results in the appearance of a nonperturbative dimensionful pa¬ 
rameter 


A = 



OO 7 

da s 

s (M) Pi&s) 



(4) 


where a s = g^/An, and (3(a s ) = da s {M)/d InM is the Gell-Mann-Low function. Furthermore, as it is well known, the 
quark mass has anomalous dimension and depends on the scale M. The renormalization -group equation for mo(M), 
the running mass, is d\nmo/d\nM = and we use the MS scheme for which (3 and 7 m are independent of the 
quark mass |27| . Upon integration the renormalization-group invariant mass is given by 

ras(M) ( \ 

m q = m oq (M) exp{ / da s } , (5) 

J p{pt 8 ) 


where the indefinite integral is evaluated at a s (M). Next we note that since free energy is renormalization-group 
invariant quantity its anomalous dimension is zero. Thus Fr has only a normal (canonical) dimension equal to 4. 
Making use of the renorm-invariance of A, one can write in the most general form 


P _ K±f( T H m n m d 

Ffi - A/( A’FT’T ) 


where / is some function. From (@)J) and (|j) one gets 


8Fr 

3(1/S 0 2 ) 


8F r dA y, 
dA d{l/gl) + 4- 


8Fr 8m q 
8m q 8(1/g%) ’ 


( 6 ) 

(7) 


8m a 


8(i/gl) 


— —Anaim, 


7 m q (a s ) 
/3(a s ) 


With the account of (|3j) the gluon condensate is given by 


<G 2 )(T, H, m u , m d ) = (4 - T-^ - 2 - ^(1 + 7m Jr 

Q 


P(a s ) 


dm. 


■)Fr- 


( 8 ) 


(9) 


It is convenient to choose such a large scale that one can take the lowest order expressions, P(a s ) 
b = (11 N c — 2Nf)/3 and 1 + > 1. Thus, we have the following equations for condensates 


(G 2 )(T, Hl m u ,m d ) = -^f(A-T-^-2H^^ 


8H 


8 

1 dm, 


■)Fr = -DFr , 


—6a 2 /27r, where 


( 10 ) 


( qq)(T,H,m u ,m d ) = 


8Fr 

8m a 


( 11 ) 


III. FREE ENERGY OF THE QCD VACUUM AT LOW TEMPERATURE IN A MAGNETIC FIELD 

The above equations enable to obtain the values of the condensates as functions of T and FI provided the free 
energy density is known. To get the latter the ChPT will be used. At low temperatures T < T C (T C is the chiral phase 
transition temperature) and for weak fields FI < n^ adr ~ (AwF ^) 2 the characteristic momenta in the vacuum loops are 
small and theory is adequately described by the low-energy effective chiral Lagrangian L e s §!• This Lagrangian 
can be represented as a series expansion over the momenta (derivatives) and quark masses 

L e g =L ( 2 ) +L ( 4 ) +L (6) + ... 


( 12 ) 
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The leading term in (|l2|) 


is similar to the Lagrangian of the non-linear sigma model in the external field 


L< 2 > 


F 2 

-^-Tr(y„U + V^U) + £i?e Tr(MU + ), 


X7 ll U = d fl U-i[U,Vj. 


(13) 


Here U is a unitary SU(2 ) matrix, F n = 93MeV is the pion decay constant, and £ has the meaning of the quark 
condensate £ = | (uu) | = | (dd) \. The external Abelian magnetic field H is aligned along the z -axis and corresponds 
to = (t 3 /2)A tl {x) with the vector-potential A M chosen as A M (x) = S^Hxi- The mass difference between the u 

and d quarks appears in the effective chiral Lagrangian only quadratically. Further, to obtain an expression for the 
quark condensate in the chiral limit we use only the first derivative with respect to the mass of one of the quarks. 
Therefore, we can neglect the mass difference between the u and d quarks and assume the mass matrix to be diagonal 
M = ml. 

At T < T c , H < p. 2 adr the QCD partition function coincides with the partition function of the effective chiral theory 

ZeS [T, H] = e -dVF eft [T,H] = Zq[h] f [DU] exp |_ f 0 dxA r d 3 xLeS ^ A ]| ( 14 ) 

J Jo Jv 

At the one-loop level it is sufficient to restrict the expansion of L e ff by the quadratic terms with respect to the pion 
field. Using the exponential parameterization of the matrix U(x) = exp{iT a n a (x) / F^} one finds 

L (2) = ^(^M 7r °) + + (<^71-+ + iA l _ l n + )(d l _ l Tr~ - iA^TT-) + M 2 7r + 7r“, (15) 

where the charged n ^ and neutral 7r° meson fields are introduced 

7T ± = (7T 1 ± in 2 )/V2 , 7T° = 7T 3 (16) 

Thus ( |l4| ) can be recasted into the form 2 

[T,H] = Z~l Z 0 [H] f [DTT°}[Dn + }[DTT-}exp{- f dx± [ d 3 xL^[n,A}} (17) 

J Jo Jv 

where partition function is normalized for the case of perturbation theory at T = 0, H = 0 

Zp.t. = [det(—<9 2 + M 2 )]" 3 / 2 . (18) 

Integration of (0) over 7r-fields leads to 

2% = Z-lZ 0 [H][det T (-dl + M 2 )]- 1 / 2 [det r (-|^| 2 + M 2 )}~\ (19) 

where — iA^ is a covariant derivative and a symbol ”T” means that the determinant is calculated at finite 

temperature T according to standard Matsubara rules. Taking (|l8|) into account and regrouping multipliers in (0) 
one gets the following expression for Z ^ F 


Z&[T,H] 


Zo[H] 


det T (—<9 2 +M 2 ) 

-1/2 

I" det(—|U M | 2 + M 2 )l 

_ det(—<9 2 + M 2 ) _ 


L det(—9 2 + M 2 ) J 


X 


rdet T (-|^| 2 + M 2 )l 
. det (—|+ Ml 


-1 


( 20 ) 


The partition function Z ^ describes charged 77^ and neutral 7r° ideal Bose gas in magnetic field. Relativistic charged Bose gas in 
magnetic field at finite temperature and density with application to Bose-Einstein condensation and Meissner effect was studied in Refs. 

@ {§,& 
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Then the effective free energy can be written in the form 


Fes(T, = In Z% = ^- 2 +F„o (T) + F„± (H) + F dia (T, H). 


Here F n o is the free energy of massive scalar boson 

d 3 p 


F, 


o{T) = T J ln(l - exp(-^p 2 + M%/T)), 


F n ± is a Schwinger result for the vacuum energy density of charged scalar particles in the magnetic field 

Hs 
/o 


Fl(m - 1 r ds r -Mls r Hs 
w ^ ^ 16 t r 2 i 0 s 3 Kmh(Hs) ^ 


L sinh (Hs) 

and Fdia is the diamagnetic free energy of relativistic charged Bose gas 


TJ'T' 00 /*oo 

F dia (T,H) = — S'] / dk ln(l - exp(-w„/T)), 
77 n =0 J ° 


( 21 ) 


( 22 ) 


(23) 


Wn = \/k 2 + M% + H(2n + 1), 

where u n are Landau levels of the mesons in constant field F[. 3 

By expanding integrand in (p2h, (p3) in the series, one obtains the following expressions 


(24) 


^71-0 — 


27T 2 


ri 

n= 1 




(25) 


and 

RT °° 00 1 / u 

Fdia = -^EV^ + H(2n + 1) X] 1*1 ( + H(2n + 1)) (26) 

71=0 fc=l ' 

where K n is the Mackdonald function. 


IV. QUARK CONDENSATE AT T ^ 0 AND H ± 0 

The free energy F^g determines the thermodynamical properties and the phase structure of the QCD vacuum state 
below the temperature of the chiral phase transition, i.e. in the phase of confinement. 

In order to get the dependence of the quark condensate upon T and U use is made of the Gell-Mann-Oakes-Renner 
relation 


F 3 M 3 = — -(m„ + rrid){uu + dd) = 2mS 


The quark condensate shift is then determined by 

AE (T,H,M n ) 


1 dF e R s 
FI dMl 


(27) 


(28) 


3 Technically, a transition for the free energy F = ^Tr In(p| + o;q(p)) from the vacuum case (H = 0, T = 0) to the case of H 7 ^ 0, T 7 ^ 0 
is straightforward. Omitting the details of the calculations, we note that, eventually, this transition reduces to the substitutions 
p 4 > ujfc = 2irkT (fc = 0, ±1,...), a* = \/p 2 + Ml - = s/p! + M% + H(2n + 1) and Tr §£ EZo EtZoo f-Z 

where the degeneracy multiplicity of H/2n has been taken into account for the Landau levels. Performing summation over Matsubara 
frequencies, we obtain (|24j). 
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Expression for quark condensate as function of T and H at M 1r ^ 0 is then given by 

(qq)(T,H,M v ) &£„o(T,M„) AS „±{H,M„) AS dia (T, H, M K ) 

(qq) S £ + S 

where the contribution of thermal 7r°-meson to the shift of (qq) is 

AS t o tm. 


A-k 2 F 2 ^ n 

11 n =1 


^ n U T 


vacuum contribution of 7r -mesons (Shwinger term) is 


AE, 


H 


dx 


(47rFV) 2 Jo 3:2 


- 1 


exp 


M 2 

_7T 


and "diamagnetic contribution” of charged Bose gas of 7r -mesons to the quark condensate is given by 

OO OO / j \ 

■EE A '« (^vm+ H (2n+i) 


ASdia H (k 


2n 2 F 2 


n= 0 k =1 
2 


In the limiting case of weak magnetic field F[ <C M 2 and low temperature T <C one gets 

AE -° = -y/8 n M : /2 r/o e~ M * IT + 0(e~ 2M ^ T ) 


and 


AS, 


(4 irF„y 


H 2 


96tt 2 F 2 M 2 


0 ( 


MF 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


Let us now consider diamagnetic term ASdia for various limiting cases. In the case of weak magnetic field, the 
interval between Landau levels ~ y/~H is much less then average thermal energy of 7r ± -mesons, y/H <C T, and we can 
use Euler-MacLaren summation formula. We then find for ASdia 


ASdia _ a AS T o 


H 2 


^8n 2 F 2 M n T 


M, 


5^ nK[ ( n ~f~ ) > ^ < T 


For the case of low temperature T <C M , one obtains 

H 2 


ASdia _ ^AS^O 


A8\/2ir 3 / 2 F 2 Mn ,2 T 1 / 2 


6 -m,/t Vli «T« M*. 


(35) 


(36) 


In the opposite limit of low temperature and strong magnetic field (at M„ ^ 0) asymptotic behavior of ASdia is 
determined by n = 0, k = 1 in the sum (B2) 


AS, 


dia 


H kJS* 


H 


HT 1 / 2 


2n 2 F 2 


T 


( 27r )3/2 F 2 (M 2 + H) I/ 4 


T <C \/~H <C M n 


(37) 


An interesting phenomenon reveals itself in the vacuum QCD phase structure under consideration. One can find 
from (|2^) such a function H(T) that the chiral condensate ( qq)(T , U, M w ) remains unchanged when the temperature 
and magnetic field change in accordance with FL* = H(T). Then Lf* is found by solving the following equation 
(qq)(T,H.,M„) - (qq) = 0 (see ©). At low temperature, T <C M w , and weak field, y/H <C M n , we have 


H*(T, M v ) = V2A 1/4 M 7 J 4 T i/4 e -M,/2T + o(e-M-IT) 

Let us now consider quark condensate as function of T and H in chiral limit AI K = 0 p2| . 


( 38 ) 
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Substituting 


into (|2^), calculating the derivative over M 2 and then taking the limit M 2 

H 


0 one gets 


ff) = <,-,)(! - | • ^5 + (4lrn)a 


In 2 — 


_H ,VH 

2tT 2 F 2 ^ T )) 


ip{\) = 


dx 


n— 0 1 


ui n (x)(exp(\u) n (x)) - 1 ) ’ 


uj n { x) = \/x 2 + 2n + l 


(39) 


Now we consider various limiting cases. In the strong field, y/ff > T (A > 1), the lowest Landau level (n = 0) gives 
the main contribution to the sum (|39|) 


<p (A » 1) 



(40) 


In the opposite limit of weak field, y/H <C T (A <C 1), the sum in J39| ) is calculated with required accuracy using 
the Euler-MacLaren formula. Furthermore, one gets the following result with the use of the asymptotic expansion of 
integral (Ehl) [fuj at A <C 1 


tp(\ < 1 ) = 


n 2 1 

YA2 


7tt 1 
24 A 


1 


C(3) 


+ -lnA + C+^A +0(A ), 


(41) 


here C = |(y — ln 47 r — |), 7 = 0.577... is Euler’s constant and £(3) = 1.202 is Riemann zeta function. Thus, one 
obtains the following limiting expressions for the quark condensate in the chiral limit in a magnetic field at T ^ 0 


(qq)(T,H) = 1 _ 1 _ + 


H 


(qq) 


3 8 F 2 (4 irF^y 


■ In 2 — 


H 3/4 T^ 2 ~i /T 

(2n) 3 / 2 F 2 


and 


(qq)(T,H) =1 _ T 2 


(qq) 


H 7sfHT 

~ — 


H 


8 F 2 (47t EV ) 2 48t tF 2 (4t tF^) 2 T 2 


Vh/t 

5 

\Th^t 

(42) 

In 

rp2 ’ 

Vh<^t 

(43) 


where A = In 2 — 8 C' ~ 4.93. 

In the framework of ChPT the quark condensate (|3!]) at H ^ 0, T ^ 0 is determined by three dimensionless 
parameters H / (AnF ^) 2 , T 2 /F 2 and A = y/H/T. The quantity A is a natural dimensionless parameter in this approach. 
The motion of a particle (massless pion) in the field is characterized by the curvature radius of it’s trajectory, 
and in the magnetic field this is the Larmor radius Rl = 1/y/H. On the other hand, there is another length 
It = 1 /T - ’’temperature length” at T ^ 0. Therefore, charged 7 r ± -mesons in magnetic field effectively acquire 
’’mass”, TO e ff = y/~H, determined by the lowest Landau level, when Larmor radius of a particle in the field is much less 
than It (A 1). Correspondingly, their contribution to the shift of the chiral condensate is suppressed by the Boltzman 
factor oc exp{— m e s/T}. In the weak field limit 7 r -mesons give standard temperature one-loop approximation ChPT 
contribution to (qq). Besides, additional temperature and magnetic corrections appear. Neutral 7 r°-meson contributes 
to ( qq)(T , H) as usual massless scalar particle. 

In the chiral limit the equation for H*(T) takes the form 


1--4A 2 ln2+^AV(A) =0 (44) 

Z7 T z 7r z 

The numerical solution of ( fi4| ) yields A* = 3.48... Thus, quark condensate stays unchanged when T and H are 
increased according to H = 12.11 • T 2 . Hence it is possible to say that the order parameter (qq) of the chiral phase 
transition is ’’frozen” by the magnetic field. 

Note that H(T c )/(4 ttF t ) 2 ~ 0.2 -C 1 at T = T c ~ 150MeV and therefore the above relations remain valid up to the 
deconfined phase transition point 4 . In the vicinity of T c the effective low energy chiral Lagrangian fails to provide 
an adequate description of the QCD vacuum thermodynamical properties, and strictly speaking becomes physically 
invalid. The following is worth noting. In deriving (p9|), at the first step the physical quantity as functions of 
where obtained, and only then the chiral limit M^ —> 0 was taken. Acting in the inversed sequence we would have 
obtained all temperature corrections to condensate identically equal to zero. This points to the fundamental difference 
of the two cases: the exactly massless particle and the particle with infinitesimal small mass. 


4 There was an arithmetical mistake in paper jiq . and the value of A* was underestimated. However, this mistake does not affect any 
physical consequences. 















V. CONCLUSION 


It has been shown in the present letter that the quark condensate is ’’frozen” by the magnetic field when both 
temperature T and magnetic field H are increased according to the obtained relation H = /7*(T, M, r ). This points 
to the fact that the direct analogy between the quark condensate in QCD and the theory of superconductivity is 
untenable. In the BCS theory the Cooper pairs condensate is extinguished by the temperature and magnetic field. The 
’’freezing” phenomenon can be understood in terms of the general Le Chatelier-Braun principle 5 .The external field 
contributes into the system an additional energy density H 2 / 2. The system tends to compensate this energy change 
and to decrease the free energy by increasing the absolute value of the quark condensate: A e v = — m\T,(H) — £(0)| < 0. 
On the other hand, if the temperature of the system is increased (by bringing some heat into it) the processes with 
heat absorbtion by damping the condensate are switched on. The interplay of these processes is at the origin of the 
above ”freezing” of £(T, H). 

In Nf = 2 QCD, the temperature phase transition restoring chiral SU(2)l x SU(2)r symmetry is a second-order 
phase transition. As the temperature is increased, the order parameter £(T) decreases monotonically, vanishing at the 
critical point T c . The existence of the jump A£(T C ) ^ 0 would indicate the occurrence of a firs-order phase transition. 
The order parameter at any non-zero magnetic held H is greater than that at H = 0; hence, £(T, H) — £(T, H = 0) > 0. 
Because of this, it may turn out that the chiral phase transition in a magnetic held becomes a hrst-order phase 
transition, and critical temperature of phase transition is larger at H 0, T C {H) > T C (H = 0). In order to clarify the 
character of the phase transition, it is necessary to investigate the behavior of the system and of the order parameter 
in the fluctuation region near T c ; however, the effects of nir interaction are substantial in this region. It follows that, 
within the approach used here, it is impossible to establish the change in the order of transition, although it seems 
highly probable 6 . This phenomenon may prove to be of use in investigating various cosmological scenarios after the 
Big Bang. For this reason, it would be interesting to continue studying the chiral phase transition in a magnetic field. 
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